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Abstract. As a first step towards a theory of differential equations involving para- 
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tions for equations of low order are given explicitly. A connection to n-generalized Fibonacci 
numbers is established. Several other classes of differential equations (systems of first or- 
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t:;^ ■ 1 Introduction 

00 ■ 

I Schwinger introduced anticommuting variables (i.e., Grassmann algebras) into physics in order 

• ' to extend his treatment to fermions of quantum fields using Green's functions and sources [1]. 

, A few year later Martin extended Feynman's path integral method of quantization to systems 

I containing fermions and needed a "classical" fermion to quantize, thereby extending Schwinger's 

■ analysis of Grassmann algebras [2]. These anticommuting variables 9i satisfying 6i9j = —OjOi - 

implying in particular that Of = Q - were used by a number of authors to develop a systema- 
^ \ tic quantization of fermions in analogy to the well-known quantization of bosons based on the 

' conventional commuting variables, see, e.g., the monograph of Berezin [3j for an early survey. 

^ Around 1970 the idea of a symmetry between bosons and fermions - called supersymmetry - 

emerged (for accounts on the early history see, e.g., [1]), providing a link between particles 
having different statistics. It took then only a few years until the commuting and anticom- 
muting variables were merged and interpreted as coordinates on a "superspace" (mainly due to 
Akulov and Volkov) and the physical superfields as "functions" on superspace (mainly due to 
Salam, Strathdee, Ferrara, Wess and Zumino). This happened in 1974 and from then on the 
idea of a physical and mathematical superworld has captured the imagination of physicists and 
mathematicians alike (see, e.g., [11 [6] for a physical discussion and [3 El [9] for a mathematical 
discussion) . In our context it is particularly interesting to note that the introduction of the anti- 
commuting (and nilpotent) variables has led to beautiful extensions of well-known mathematical 
structures (supermanifolds, supergroups, etc.). 

In a rather different line of thought Green introduced in 1953 parastatistics into quantum 
field theory by considering certain trilinear commutation relations, generalizing the quantization 
schemes for bosons and fermions [10]. With hindsight it is not surprising what the next steps 
should have been in analogy to Fermi-Dirac statistics (i.e., fermions): The introduction of the 



2 



T. Mansour and M. Schork 



"classical paraobjects" which yield upon quantization the operators satisfying the trilinear com- 
mutation relations (and are nilpotent in higher order). This should motivate the introduction 
of para-Grassmann algebras and then, in close analogy to above, a merging of commuting and 
para-Grassmann variables as coordinates on a parasuperspace and describing physical parasu- 
perf ields as functions on this "space" . From a more mathematical point of view one would hope 
that these parastructures would again yield fascinating new extensions of well-known structures. 
In reality the development took some time; although there had been precursors [21 [11] the real 
development took off around 1979 with the work of Kamefuchi and coworkes [121 [T3| I14j . In 
their work the focus was on the quantum field theory of particles satisfying parastatistics, see, 
e.g., the early monograph [15j . Para-Grassmann variables - satisfying [6j, Of^]] = as well as 
[6i-i^,6i2, . . . , di^]+ = where m > p + 1 and [.,...,.]+ denotes the completely symmetrized pro- 
duct (note that this implies 9f~^^ = 0) - were mainly considered as a convenient tool, but already 
in [13] parasuperfields depending on conventional and para-Grassmann variables were introduced 
in close analogy to superf ields. One of the first applications were in a model of one-dimensional 
supergravity interacting with "matter" represented by para-Grassmann-variables [16] and in 
a para-Grassmann extension of the Neveu-Schwarz-Ramond algebra [T7|. Further applications 
have been in a parasupersymmetric extension of quantum mechanics [18^ [T9| [20l [211 I22j . in 
particular in generalizations of the Pauli and Dirac equations [231 121] to describe particles with 
spin s > ^, as well as in the closely related fractional supersymmetry [25l [26l [27l [28l [29] . 
Low dimensional field theories, in particular conformal field theories, have also been discussed 
[30l[3ll[32l[33l[34l[35l[36l[37l[38l[39] and a generalization of the Virasoro algebra has been discove- 
red [401 [¥T] . In particular in this approach the geometrical point of view of merging the ordinary 
variables with the para-Grassmann variables into a set of coordinates on a parasuperspace and 
considering the fractional supersymmetric fields as functions on this space has been used in close 
analogy to the well-known superspace and superfields. From a different point of view a connec- 
tion to deformed algebras or quantum groups was reported [42l HH HH |45l HB] |47l US] [49] . Of 
course, para-Grassmann variables have also been considered in their own right [501 1511 [521 153] . 
Many more references can be found in [54] which surveys also some of these developments. 

Let us point out that the para-Grassmann variables considered originally (satisfying trilinear 
relations) lead quickly to rather messy structures and calculations. Therefore, a generalization 
of the Grassmann variables was introduced by Filippov, Isaev and Kurdikov in [411 [42l [51] 
satisfying 9^ = as well as certain bilinear commutation relations of the form 9idi = ep+i 9j9i 
(for i < j). In most of the more recent works (i.e., after 1992) these generalized Grassmann 
variables - also called para-Grassmann variables - are used. Of course, in the single-variable 
case these different approaches reduce to the same sole condition 0^+^ = and the correspon- 
ding structure including a parasupercovariant derivative has been introduced in [30] (a formal 
variable A satisfying A^'"'"^ = was already considered by Martin in [2]). 

From a more mathematical point of view one may hope that there exists an associated Zp+i- 
graded mathematics, yielding for p = 1 "ordinary" supermathematics. Some results in this 
direction have already been obtained in the literature cited above. In particular, in the paper [53] 
by one of the authors many basic elements of para-Grassmann calculus (involving a single 
para-Grassmann variable) are collected and some simple observations concerning differential 
equations involving a para-Grassmann variable were obtained. It is the aim of the present paper 
to start a systematic treatment of differential equations involving para-Grassmann variables. 

Before turning to a more detailed description of the present paper we would like to point out 
that in [38] a differential equation containing a para-Grassmann variable is solved in connection 
with deformed coherent states. To the best of our knowledge this is the only paper where such 
a differential equation is solved explicitly (more details can be found below). 

Now, let us describe the contents of the present paper in more detail. In Section [2] the 
definitions and basic properties of para-Grassmann calculus are recalled and the notation is 
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fixed. In Section [3] some of the simplest differential equations are solved and compared to the 
usual "bosonic" differential equations. The linear differential equations of second order (with 
constant coefficients) are considered in Section [H In the direct approach chosen a recurrence 
relation has to be solved which is closely related to Fibonacci numbers. The results obtained 
are extended to linear differential equations of arbitrary order in Section [5j However, since the 
method used for the second order case turns out to be rather cumbersome in the general case 
(involving recurrence relations related to n-generalized Fibonacci numbers) we follow a slightly 
different approach. The main difficulty consists in getting control over the case where degen- 
eracies appear in the associated characteristic equation. In Section [6] the general system of first 
order is considered. It is rather simple to describe the solution of this system abstractly but 
when considering concrete examples tedious calculations result, too. The peculiarities of the 
linear differential equation with variable coefficients are sketched in Section [7] by considering an 
example. In Section [8] a very simple nonlinear differential equation is treated. It turns out that - 
in contrast to the linear differential equations where a close analogy to the "bosonic" case holds 
true - due to the nilpotency of the para-Grassmann variable the equations show a behaviour 
reminding one also of "bosonic" partial differential equations. Finally, some conclusions are 
presented in Section [H 



2 Some formulae of para-Grassmann calculus 

In this section we briefly recall some basic facts about para-Grassmann calculus following [5^ 
(see also |30 1 l42 l H H l51]). For a given natural number p > 1 the nilpotent "coordinate" vanishes 
in {p + l)-th order, i.e., 6'^"'"^ = (with 0' / for < / < p); the case p = 1 reproduces the well 
known Grassmann variable appearing in super symmetric theories. Due to our assumption there 
will be the linearly independent (over C) elements {1,9,9'^, . . . ,9^}. In analogy to "ordinary" 
supermathematics we introduce "parasuperfunctions" /(x, 9), which can be expanded as a finite 
series in 9: 



fix, 9) = foix) + Mx)9 + • • • + U{x)9P. (2.1) 

The functions /fc(x) will be called the components (or sometimes coefficients) of f{x,9). We 
introduce a derivative d = dg satisfying the following commutation relation: 

d9 = q9d, q = e^. (2.2) 
In the following we use the standard notation 

[n], = l + q + ...+ ^ 1^, [n],\ = [n]Jn - 1], • • • [1], 

with [0]g = [0]g! = 1. It is easy to check that the commutation relation ()2.2p implies 89 = 1, 
89"^ = (1 + q)9, and in general 

Note that d^^^ = 0. For parasuperfunctions we have (at least) two derivatives, the "normal" 
derivative dx = and d acting on the nilpotent variables. We introduce a parasupercovariant 
derivative by [30l SH El EI] 



(2.3) 
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It is straightforward to check that for / as in (|2.ip one has (for < r < p) 

^'/(.. ») = 1 |^JL_a,A,.„.-(.-..« + ± ^A(^)«-'. (2.4) 

Note in particular that the parasupercovariant derivative X> is a (p+ l)-th root of dx in the sense 
that [301 da mi EH 

vp+'f{x,e) = dxf{x,e). (2.5) 

Due to this relation there exists a certain "periodicity" in the derivatives "D^f of a function 
/ where now s is an arbitrary integer. Since s can be decomposed in a unique way as s = 
s'(p + 1) + s" with < s" < p + 1, one may use (p^ and (p3]) to obtain 



+ E|F^8J'A(.)e'-"". (2.6) 



It is also possible to introduce a formal "antiderivative" by letting 

Vfix, 9) := \p],\ f fp(t) dt + Ux)e + -^Jii^W^ + ■■■+ [^/p-i(^)^'' 

so that TD^^Vf = /; in particular, = 9. More generally, iterating this allows one to 

introduce an operator Let us also introduce an exponential function in a bosonic and 

a para-Grassmann variable by setting 



an 

Then it is an easy calculation to verify that it has the defining property of the exponential 
function, i.e., 'Deq{x;9) = eq{x;9). Under a scaling of the arguments we obtain the following 
nice expression 



Veq{pP^^x;(i9) = l3eq{l3P+^x; (59). 
This implies VP+'^eq{(jP+'^x; (56) = dxeg{(3P+^x; (36), as it should according to ([23]) . 

3 Some very simple differential equations 

In this brief section we consider some simple differential equations. Most of these results can 
be found already in [53] but for the convenience of the reader and the systematic exposition we 
have included the results here. Now, let us consider the very simple differential equation 

V'^fix,6) = 0. (3.1) 

Writing s = s'(p+l) + s" with < s" < p+l, we may use (j2.6p to obtain the following conditions 
on the components fk{x): 

di+^h{x) = 0, < A: < s" - 1, dih{x) = 0,s"<k< p. 
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Thus, the first s" components are polynomials of degree s', whereas the remaining components 
are polynomials of degree s'—l. In particular, there are s"{s'+l)+{p+l—s")s' = s'(p+l)+s" = s 
free parameters for a solution / of (j3.ip . showing that the kernel of has dimension s. In the 
case s = 1 this means that the only solution of Df = is given by f{x,9) = c for some c E C. 
Thus, the "parasupercovariant constant function" (meaning that Vf = 0) is indeed given by 
a constant value f{x, 9) = c. In the case s = 2 the "Laplace equation" P^/ = in one dimension 
has two linearly independent solutions. Consider the inhomogeneous version of (j3.ip . i.e.. 



V'f{x,e)=g{x,e). (3.2) 

A particular solution is given by 'D~^g. For any solution /hom of the homogeneous equation the 
sum /horn + 1^~^9 is a solution of the inhomogeneous equation. Thus, the set of solutions is 
an affine space (as in the usual "bosonic" case). Let us summarize these observations in the 
following proposition. 

Proposition 3.1. Let the differential equation (13. ip with s = s'{p + l) + s" where < s" < p+1 
be given. Any solution of (13. ip can be written as 




fhom{x, 9)= y ^ {y ^ Ck^rX" )6 + / , S Ck^rX 



where Ck^r S C. Thus, the space of solutions is a complex linear space of dimension s. The 
associated inhomogeneous differential equation (j3.2p has the particular solution 'D~^g. For any 
solution fhom of the homogeneous equation the sum /hom + T^"9 "is a solution of the inhomoge- 
neous equation. Thus, the set of solutions of (j3.2p is a complex affine space of dimension s. 

Now, let us turn to the "eigenvalue problem" 

vf{x,e) = xf{x,e) (3.3) 

for A € C. From above we know that f{x, 6) = Ceq{XP~^^x; X6) is a solution. To find out whether 
there exists another solution, we write f{x,d) = 'Yl\=o fk{x)0^ . Inserting this into (13. 3p gives 
for A; = 0, ... ,p — 1 the recursion relation fk+i{x) = j ^^^j^ i fk{x) as well as dxfo{x) = \[p]q\fp{x). 

Using the recursion relation repeatedly yields fp = i^pj/o- Using now the second equation yields 
the differential equation dxfo{x) = X^~^^fo{x) for /o with the solution fo{x) = Ce^^^^^ where 
C G C is a constant. From the recursion relation it follows that fk{x) = '^\iCe'^^^^^ and, 

consequently, that f{x,6) = Ceq{XP^^x; X6) is indeed the general solution of the differential 
equation (|3.3p . Using (|2.6p . it is straightforward to consider the slightly more difficult eigenvalue 
problem 

V'^fix,e) = Xfix,e), sGN. (3.4) 
Namely, let A~ be one of the s roots satisfying (A~)^ = A. All s roots are then given by 

{x-s,nx^,...,n'-^x-^}, 

where {I, fi, fi"^ , . . . is the cyclic group of order s (isomorphic to Zg) consisting of the s 

roots of unity. Since each function eq{{^^ X^Y'^^x; {fi^Xs)6) solves (|3.4p . the solution of (j3.4p is 
given by f{x,9) = X^^^q C'fceg((/i'^As )p+^x; {fi''X'i)6). The space of solutions has dimension s, 
as expected. Let us summarize these observations in the following proposition. 
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Proposition 3.2. Let the differential equation (|3.4p be given. Let he one of the s roots 
satisfying {\~y = A and let {1, ^u, /i^, . . . , ^u'*""'^} he the cyclic group of order s. Any solution 
of (j3.4|) can he written as 

/(x,e) = X;C,e,((/A^r+^x; (/A^)0), 

where G C. In particular, the set of solutions of ()3.4p is a complex linear space of dimension s. 

Up to now, we have only considered equations wliere a simple ansatz using the exponential 
function ()2.7p yielded all solutions. In the following sections more complicated expressions in T> 
will be discussed. Due to the lack of a general product rule these cases will turn out to be more 
difficult and it is unclear what a "good" ansatz should be. This will appear already in the case 
of linear differential equations of second order discussed in the next section. 

Remark 3.1. Let us denote the set of parasuperfunctions / = f{x,9) by A (where we assume 
for ease of presentation that the components fk are "sufficiently nice", e.g., in C°°). Clearly, it 
is a linear space (over C) and if we also consider as product of elements f,g & A the pointwise 
product (i.e., {fg){x, 9) := f{x, 6)g{x, 6)) it is even a commutive ring with unit (having nilpotent 
elements). The parasupercovariant "derivative" T> : A ^ A is, however, not a derivation in this 
ring, i.e., 'D{fg) ^ T>{f)g + fT>{g). Thus, the structure {A, is not a differential ring as con- 
sidered, e.g., in [55] and the literature given therein. It is tempting to introduce a new product * 
such that {A, *,T>) becomes a differential ring, i.e., P(/ * g) = 1^{f) * g + f * T^{g)- However, 
these conditions on the new product seem to be rather involved and a natural interpretation is 
lacking. 



4 The linear differential equation of second order 

In this section we will discuss the general linear differential equation of second order. More 
precisely, we will show the following theorem: 

Theorem 4.1. Let ci,C2 E C; the general solution of 

{V^ + ciV + C2)f{x,e) = (4.1) 



may be obtained as follows. Define A-t := — y ± ^ y^cf — 4c2 . In the non- degenerated case where 
A+ 7^ A- the general solution of (|4.ip is given hy 

fix, 6) = CiCqiX^^^x; X+e) + C2eg(A!+^x; X.9) 

with arbitrary Ci € C. In the degenerated case where cf = 4c2 we abbreviate a := — y [thus 
C2 = oP'); the general solution is in this case given by 



fix, 6) = Cieg{aP+^x; a0) + C2 (p + l)aP xeq{a^+^ x; a0) + e""+'" Yl lu ■ (^.2) 

V k=i ^''J^- 



Note that making an ansatz of exponential form (as we did in the last section) yields in the 
non-degenerated case both solutions and in the degenerated case the first solution. Before giving 
a proof of the theorem, we single out the cases of smallest p explicitly. 
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Corollary 4.1. Consider the degenerated case (T)"^ — 2aD + a^)f{x,9) = 0. ^ first solution is 
always {i.e., for arbitrary p) given by fi{x,9) = eq{a^~^^x;a9). A linearly independent solution 
is given in the case p = 1 {i.e., for an "ordinary" Grassmann variable) by 

f2{x, 6) = 2axe"'^ + (l + 2a^x)e''''''e. 

In the case p = 2 a second solution is given by 

/2(rE, 0) = 3a\e"'^ + (1 + 3a3x)e°'^e + ^ (2 + 3a=^x)e"'^02^ 

Let us now turn to the proof of the theorem. Instead of just inserting the claimed solutions 
and see whether they fulfill the differential equation we will give a longer proof and show how 
these solutions are found. 

Prof of Theorem Sm Let f{x,e) = fo{x) + fi{x)e + ■■■ + fp{x)9P be a solution of KTh . 
Inserting this and comparing coefficients of equal power in yields the following system of 
equations: 

^'^'^'^ = -'''JkTWj''-'^'^ - [fc + 2Ufc + 1]/ ^' 0<k<p-2, (4.3) 

Oxfo{x) = -Cl\p]qlfp{x) -C2\p- l]q\fp-l{x), (4.4) 

dxfi{x) = -cidxfo{x) - C2\p]q\fp{x). (4.5) 



Note that in the case p = 1 the recursion relations (j4.3p are vacuous. Only the differential 
equations (I4.4p - (l4.5p remain and can be cast into the form 



dxfi{x) J V C1C2 (C1-C2) J V fi{x) J 

Let us return to the case of arbitrary p. Here we will first solve the recursion relations until we 
arrive at a similar system of differential equations for the components fo{x) and fi{x) which 
can be solved. To solve the recursion relations, we introduce for k = 0, . . . ,p — l the vectors and 
matrices 

Vk{x) := ( /'^f. ) , Tfe+i := ( % \ ] . (4.7) 

V Jk+l[x) ) y -lfc+2|,lfc+l|, -IFf2i; j 

This enables us to write the recursion relation (j4.3p in the compact form v\^j^\{x) = Tk+iVk{x). 
Iterating this yields 

Vk+i{x) = Tk+iTk ■ ■ ■ TiVo{x) = Fk+iVo{x). (4.8) 

Since Vp-i{x) = Fp-iVo{x), this allows us to express fp-i and fp through /o and /i once we 
have determined Fp_i. Inserting this in (j4.4p and (j4.5p will give a system of ordinary differential 
equations for /o, /i. We will now determine for all k since after having determined /o and /i 
the remaining components fk are calculated using ()4.8p . Let us write the matrix Ffe explicitly as 



/ fW j-(2) 

fW r{2) 
Jk Jk 

Note that we can then write 



(4.9) 



AW = ||/o(.) + ("0) 
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The relation Fk+i = T^+iFfc is equivalent to 

/S2 = -ci/Si-C2/f, i = l,2. (4.11) 

In general, a solution of a generalized Fibonacci sequence fk+2 = —cifk+i — C2fk is obtained as 
follows: The generating function + cix + C2 has zeroes A-t = — ^ i \ \/c\ — 4c2 , so that by 
Binet's formula the general term is given by = aiA^ + 02^^, where the coefficients a\, 02 are 
determined using the (given) initial values /o, f\. Note that we have the following relations for 
the coefficients q: 

ci = -(A+ + A_), C2 = A+A_. (4.12) 

Let us now assume that A_|_ = A„, i.e., = 4c2. We define ci = —2a, so that C2 = a^. In 
this case we obtain (for k = 0, . . . ,p — 1) 



iFfiy [k+i]qi 



(4.13) 



Recalling ()4.8p . we may now use this for /c = p — 1 to express fp, fp-i through /o, /i- Inserting 
this into (14. 4p and (14. 5p gives the differential equation ^0(2;) = Avq{x), where the matrix A is 
given by 



A 



-paP+^ {p + l)aP 
-(p+l)aP+2 (p + 2)aP+i 



A has only one eigenvalue A = aP~^^ with eigenvector (1 a)*. Changing to the appropriate basis 
allows us to write 

± ( foix) \ / aP+^ ip+l)aP \ ( /o(x) 
dx\ h{x) + aUx) ) V aP+i ) \ h{x) + afo{x) 

which implies that there exist Ci , C2 S C such that 

/o(x) = (C7l+x(p+l)aPC2)e"'+'^ 
/i(x) = (aCi + x(p + l)aP+^C2 + C2)e"""''". 

On the other hand, by (j4.13p we obtain that 



h(y^ 1 



for A; = 2, 3, ... ,p. Hence, 

f{x,e) = fo{x) + fi{x)9 + }_^h{x)e'^ = fo{x)2^^jj/^ + ^^^ — -^22 



[k]g\ a ^ [k] 

Cieg{aP+^x; aO) + C2 ( (p + l)aP xe^ia^-^^ x; aO) + e"'^'" V ^" ) , 

V k=i ^''J'?- / 



which completes the proof of this case. 
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We now assume that A+ 7^ A_. Thus, we obtain the two generalized Fibonacci sequences 
k^)km and (/f 



(/fc^^)fceN and (/P'')fceN which satisfy the same recursion relation (j4.1ip but with different initial 



conditions. More precisely, one has for the first sequence /^^^ = and /^^^ = -C2, while 

(2) 12) 

the initial values for the second sequence are given by /g =1 and = — ci. Defining 
7 := \/c1 — 4c2 , we can write the general terms of these sequences with the help of Binet's 
formula as 

h - — /^fc = — (A_ - A+j, 

/f4^. + |^«.4(l-^)Ai4(l + ^)A'-. (4.14) 

thereby defining implicitly //fc, Uk- In particular, using this for k = p — 2 in (j4.8p yields 

.(1) .(2) ,(1) .(2) 

Combining (j4.4p . (j4.5p and (j4.15p gives the differential equation 
where the matrix i? is given by 

[(4- C2)4L\ + C,C24i), (4 - 02)4% + ^1-24^^2 J ■ 

Note that this is the generalization of (|4.6p for arbitrary p. Using (j4.1ip as well as (j4.14p . we 
obtain that 

B = ( \ =- ( ^"^^^ -/^p+i 

y ^^P+l 2^P+l + I^^P+l ) 7 V C2l-lp+l -fJ-p+2 

In order to find the eigenvalues of the matrix B let us first find the characteristic polynomial 
of B: 

det{B - xl) = x^ - ^{c2fip - fJ.p+2)x + ^ (/Xp+i - /.ip/ip+2) • 

Recalling (j4.12p and 7 = A+ — A_, we obtain that 

det(5 - xl) = x"- vp+ix + c^+^ = {x- A^+') {x - A^+^) . 

Hence, the eigenvalues of the matrix B are A^^^ and A^^^ with the eigenvectors (1 A_)* and 
(1 A+)*, respectively. Therefore, there exist constants Cj G C, i = 1,2, such that 

In the case p = 1 we are finished since all components fQ{x) and have been determined. 

For arbitrary p we use ()4.10p and (j4.14p to get that 
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for all /c = 2, 3, . . . ,p. This in turn yields 

p 



fix, e) = Y^ h{x)e^ = h{x) + fi{x)e + J] -^(c2^,_i/o(x) - 

fc=0 fc=2 ^ 



k=2 '"'I- ' k=2 L'^J?- 



A+/o(x)-/i(x)^(A„0)'= , /i(x)-A_/o(x),^(AH.e)^ 

t;(7 -r 

Using (I4.16|) . we find that 



/o(^) + = C7ie^-"''^(1 + A_0) + C2e^+"''"(1 + A+0) 

as well as A+/o(3;) — fi{x) = jCie - ^ and fi{x) — A„/o(x) = 7026 + ^, implying 
fix, 9) = Cie^-^'"(1 + X.e) + C72e^+^'"(1 + X+O) 

k=2 ^ k=2 ^ 



fc=0 ' k=0 



= Cie5(Ar X, A_0) + CiSqiX^^'x, X+O), 
which completes the proof. 



5 The linear differential equation of arbitrary order 

In this section we solve the linear differential equation of order n, where we assume for ease of 
presentation that n < p (for n > p + 1 one has to consider (12. 6p instead of (12. 4p ). Thus, we are 
considering the equation 

iV" + ClP'^-^ + • • • + Cn-lV + Cn)fix, 6) = 0. (5.1) 

Observe that the restriction n < p means that for p = 1 we only consider n < 1 - but this case 
was already considered in the above sections. Thus, the case p = 1 doesn't lead to new cases here. 
At the end of the section we consider briefly what happens when n > p (see Proposition 15. ip . 

The strategy for the solution of the above differential equation consists of two steps: In the 
first step we consider the linear differential equation associated to one degenerated eigenvalue 
and in the second step we assemble these results for the differential equation (jS.ip . The second 
step is only of a formal nature and the main work consists in finding the solutions of the 
degenerated differential equation. 

In the following we will need multi-indices. Therefore, let us fix the notation. We write 

rris := (mo, mi, . . . ,ms) 

where niQ > and m^ > 1 for i = 1, 2, . . . , s are natural numbers. We also introduce the length 
of rris by 

s 

\ms\ := y^mj. 
i=0 
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Lemma 5.1. Let f{x,d) = Y1\=q fk{x)(^^ be a solution of the differential equation 

(P- A)"/(x,0) = 0, l<n<p. 
Then there exist Ci , . . . , C„ € C such that 



where m^ = (?no, mi, . . . , rus) with rriQ > and mi > 1 for all i = 1, 2, . . . , s. 

Proof. We proceed the proof by induction on n > 1. For n = 1, the lemma gives ~ due to 
s = and, hence, mo = - that fk{x) = C'le'^^'^^^j^ which agrees with the solution of (j3.3p . 
Thus, the lemma holds for n = 1. Now, assume that the lemma holds for n, and let us prove it 
for n + 1. Let us write 

{V - Xr+'f{x,9) = {V - \r9{x,e) = 

with g{x, 6) = {T> — X)f{x, 6). Using the induction hypothesis (for g{x, 9)), we obtain that there 
exist Ci , . . . , C„ G C such that 



I — n I — n L iq „_n I i„ i \ '-'/ ,_i 



On the other hand, recalling g{x,9) = (T> — X)f{x,9) and using (j2.4p . we obtain 

fc=o ^^J"?- fc=i -^J''- fc=0 

By comparing the coefficients of 9^ of the above two equations, we obtain that 



mi / \ s\ 



9p{x) = -^^dxfo{x) - Xfp{x), 



gk{x) = [k + l]Jk+i{x)-Xfk{x), l<k<p-l. (5.3) 
Therefore, 

p-i p 
d.foix) - X^-^'foix) = d.foix) + X^X'\p-i]g\fp.i{x) - aE A> - Vfp-^{x) 

i=Q i=0 

P 

Noting dx{e~^''^^^ fo{x)) = e~'^''^^^{dxfo{x) - A^+-^/o(x)), we conclude that 
dxie-^'^'^foix)) = e-^^^'-EA>-i]g!5p-.(x). 

i=0 

Inserting the expression for gp^i{x) given above, this implies 

e'^^'^dAe-'^^'^foix)) 
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xp- 



n-l 



i=0 



s=0 I |ms|<n-l 



mo 



p n— 1 



1=0 s=0 I |ms|<n-l 



mo 



\ 1 1 tin 

p + l\ [ (APx)^ 



m, / si 



Thus, 



p n—l 
1=0 s=0 



^ |ma|<n— 1 



p — I 

mo 



n 



m. 



Using 



(A^x)^ _ 1 ({XPxf^^ 



s! AP V (s + 1)! 
we have found 

p n— 1 
i=0 s=0 

or, 

n— 1 f p 



C'n-|m,|A'' 
^ lms|<n-l 



p — I 

mo 



n 



p + A I (APx)^+^ 
j f (s + 1)! ' 



e-^^^^^/o(x) = 6 + j; 5; C„„|^^|A^ 

s=0 I i=0 |ms|<n-l 



p — i\-A-/p+l\ I (A^x 



mo 



n 



ma 



\s+l 



(s + 1)! 



for some constant 6 G C. Thus, using (m*) ~ (m"+i)' ^^^^ 



'^/ow-»+Ej E ^.H.,-.,A-"-'«'(„;Vi)ri('',:; 

s=l I |m,_i[<n-l V U ' / j=l \ J 



(A^x)^ 



Introducing a new multi-index = (no, ni, . . . , Ug) by setting 
no := 0, ni := mo + 1, := mj_i, 2 < i < s, 

we find |ms_i| = |ns| — 1 and, consequently. 



/o(x)=6+^ c„_i„^i+iA^-i"^i n 

s=l I |ns|<n;no=0 j=l 



p+l\ \ {XPx 



It is clear that we can consider the constant 6 as a summand corresponding to s = 0. Switching 
in notation back from to uig, this yields finally the result 



s=0 I |ms|<n;mo=0 j=l 



p+l\ I (APx) 



mi 



This is indeed the form asserted in (|5.2|) for /o(x) (note that in the formula there one has - due 
to = - that mo = holds true). 
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Now, let us find fkix) with k > 1. By an induction on k using ()5.3p we get that 

Substituting the expression of gi{x) shows that 

-A^+ix..N_ A'' ^^ J _, , -A- /p + 1\ I (A^x) 



e 



\ K I / I 1 \ 

■Aw=^E E f^»«-in..iA-'-'n(0 

^ J"?' s=0 l\ms\<n;mo=0 j=l ^ ^ ^ 

*Swl{,j_.-~'--'"'(.:.)n(-:;))'^ 



E E ^...H».A-'(:jnc:;) 

s=0 l_|m,|<n;mo=0 ^ j=l ^ ^ 



s\ 



It follows that 



which completes the proof. 



Using the expression given in (|5.2|) for fkix) (see Lemma l5.ip together with f{x,9) 
Yl^k=o fk{x)0^ , we obtain the following result. 

Theorem 5.1. Let f{x,9) = Y^^k=Q fk{x)0'^ be a solution of the differential equation 

(2?-A)"/(x,^) = 0, l<n<p. 
Then there exist Ci , . . . , C„ G C such that 



f{x,9) = e' 



k=0 {s=0 \|ms|<n-l ^ j=l 



m J j s\ I [/c]g 



where nig = (mo, mi, . . . , rris) with rriQ > and > 1, /or all i = 1,2, ... ,s. 

Example 5.1. Let us consider the case n = 2 (the case n = 1 was already considered above 
during the induction) . The formula given in Theorem 15.11 yields for s the two summands corre- 
sponding to s = and s = 1. Let us consider first s = 0. It follows that |mo| < 1 and, therefore, 
that mo = or mo = 1. The entire inner sum yields thus C2 + Ci\~^k. In the case s = I the 
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only possibility is mi = {mo, mi) = (0,1) with |mi| = 1. The inner sum yields in this case 
Ci{p + l)XPx. Thus, 



/(x, e) = e^"^'" Yl {^2 + CiX-'k + Ci{p + 1)A%} 



k=0 



{p + l)XPxe,{\P+^x; X9) + e^'^'" , j . 

This is exactly equation (j4.2p given in Theorem 14.11 for the case n = 2 (with degenerated 
eigenvalue) . 

Using the linearity of the general case, equation (jS.ip and Theorem 15.11 we get the following 
theorem. 

Theorem 5.2. Let the differential equation (|5.ip with 1 < n < p be given. Assume that the 
different roots Ai, . . . , Am, of the characteristic polynomial x" + Cix"~^ + • • • + Cn-ix + Cn have 
multiplicities ni, . . . , {with ni + • • • + = n and nj > 1). Then the solution f{x, 6) of the 
differential equation (jS.ip can be written as 

(i) ^s-~\m4 k\'fp+l\\{X^xr\{X. 



i=l k=0 I s=0 \|m.,l<n,-l '^'^^1=1^ - I I L jg 

where c|*\ . . . , Cn-"* S C are some constants {there are n of them) and m^ = (mo, mi, . . . , m, 
mo > and irii > 1, for all i = 1, 2, . . . , s. 

This theorem immediately implies the following corollary. 



Corollary 5.1. Let the differential equation (|5.ip with 1 < n < p be given. The set of its 
solutions is a complex linear space of dimension n. The set of solutions of the corresponding 
inhomogeneous differential equation 

{pn + ciP"-l + • ■ ■ + Cn-lV + Cn)f{x, 0) = g{x, 9) 

is an n-dimensional complex affine space. 

Proof. It only remains to check the last assertion. The argument is the same as in the usual 
case: If fi{x,9) is a particular solution of the inhomogeneous differential equation then any 
function of the form fhom{x,9) + fi{x,9) - where fhom{x,9) is a solution of the homogeneous 
differential equation - satisfies the inhomogeneous differential equation. ■ 

Remark 5.1. An approach similar to the case n = 2 is also possible in the more general case 
where 2 < n < p. Inserting f{x, 9) = fo{x) + - ■ ■+fp{x)9P into (jS.ip yields after a straightforward, 
but slightly tedious computation in close analogy to ()4.3p - ()4.5p the following system 

/fc+n(x) = -2_]ci J^' fk+n-i{x), 0<k<p-n, (5.4) 

j^i r + "-Jg- 

» f f \ - Sr^ [k + n-p-l-i\q\ 

OxJk+n-p~l\X) — — y Ci — - — ■ -T— — OxJk+n-p~l~i\X) 

^-^ [fe + n— p— IJq! 



i=\ 
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E < 

=k+n—p 



[k + n 



[k + n — p- l]q\ 



fk+n—ii.-^\ 



p — n + l<k<p. 



(5.5) 



Note that this system constitutes indeed p + 1 conditions as it should and that this system 
reduces for n = 2 to the one given in (j4.3p - (j4.5p (here the recursion relations (j5.4p are also 
vacuous in the case p = 1 and only the system of differential equations remains - as in the 



case n = 2). Considering the recursion relation and abbreviating ^i{k) :- 
write (15.41) as 



[k+n- 



-, we may 



/ 



fk+i{x) 



fk+n{x) 



1 





1 



1 





\ ln{k) 7n-l(fc) 



\ 





1 



( fk{x) 



\ fk+n-lix) 



Let us denote the above n x n matrix by T^'^j^ and the vector on the right hand side by 

(rf\ _____ 

Vk{x); observe that T^t|_i equals T^+i from (j4.7p . We may, therefore, write the above equa- 



.(2) 



tion as Vk+i{x) = T^'^^t;fc(x). 



This implies that Vp-n{x) 



- p—n 



Let F 



(n) 
fc+1 



(n) 



As in the case n = 2 considered in the last section, we may reformulate the re- 



cursion F 



(n) 
fc+1 



-,(n) TT,(n) 



-k+1^ k 



F, by considering appropriately normalized entries of the matrices F 



(n) 



More precisely, we introduce in analogy to ()4.9p the entries fm by 



J k—n+i 



" y[k-n + i + l]q\^ 
This implies for each j with j = 



1,... 



n the recursion relation 



fU) 

J k+n 



(5.6) 



Thus, for each j there will be an n-generalized Fibonacci sequence {fj^^^)keN satisfying (|5.6p [56l 
[57] : they are differing only in the initial values (/o"'\ . . . , fn-i)- The initial values of {fl''^)keNi 
are exactly the normalized entries in the j-th column of T^^^ . This is the generalization of (j4.1ip . 
To proceed further in close analogy to the case n = 2 one denotes the zeroes of the generating 
function + cix"~^ + ••• + £„ by Ai, . . . , A„, where some of these zeroes may occur with 
multiplicity greater than one. By using the appropriate generalization of Binet's formula (cf. [561 
ETlllH]) one can write the general solution as fk = aiA^ + - ■ ■ + anA^, where the are determined 
by the initial values. However, the resulting expressions become rather unpleasant so that we 
won't consider this procedure further. It is, however, clear that this will lead in the second step - 
after inserting the resulting relations into (j5.5p - to a system of linear differential equations which 
has then to be solved. 



In the above considerations concerning (j5.ip we have always assumed that 2 < n < p. Recall 
that the motivation for the restriction n < p comes from the periodicity ()2.5p which shows that 
in this case derivatives of at most first order in the components /fc(x) appear, see ()2.4p . whereas 
in the case n > p + 1 derivatives of higher order appear, see (j2.6p . In the remaining part of the 
section we consider the case where an arbitrary n is allowed. It is to be expected that (15. ip is 
equivalent to a system of ordinary differential equations of higher order. 
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Proposition 5.1. Let the differential equation ()5.ip be given where n = n'[p + 1) + n" with 
< n" < p+1. This differential equation for f (x , 9) = fo{x)+ fi{x)9 + - ■ ■ + fp{x)9P is equivalent 
to the following system of differential equations for the components fk{x) 

n" n' — l p 

^Cn"-u[l'W-d^' fuix) = - ^ ^C(^n'-i,){p+l)+{n"-u)Mg-9xfu{x), 
u=0 fi=0 u=0 

n" n' — l ^ p-l 

^Cn"-u[l + l^]qld^' fl+u{x) = i ^ C(^n' -t,){p+l) + {n" + I'W'd^ fl+Ax) 

u=0 11=0 ^ u=0 



+ Yl (^{n'-ii){p+i)Hn''->^)[i±Ji]g-9xft^{x)l 

,,— „ ; I 1 / 



u=p-l+l 

n" p-l 

J2 C^"-u\l±J^Wd:'+'fl+^ix) = -Yl ^n"-u[l + l^]g\d:' fi+,{x) 
;/=p+l-Z u=0 

n'-l / p-l 

- X] I y^.^(n'-fJ.)(p+l)+(n"-iy)[^ + t^W-dx fl+ujx) 
11=0 \ u=0 

P \ 

^{n' — ^){p+l)+(n" —u) 
v=p-l+l J 

where we have introduced the abbreviation x := x—{p+l) as well as cq := 1. The second equation 
holds for all I with 1 < I < p — n" and the third equation for all I with p — n" + 1 < I < p. 

Before turning to the proof we would like to point out a few facts. Note first that there 
are in total p + 1 equations which have to be satisfied, as it should be. As a second point 
note that there are n" differential equations of order n' + 1 as well as p — n" + 1 differential 
equations of order n' . The case n < p considered before corresponds to n' = and n" = n: The 
above proposition yields that there are n differential equations of order 1 as well as p — n + 1 
differential equations of order - which is indeed the case, see ()5.4p - ()5.5p . On the other hand, 
if n = n'{p + 1), i.e., n" = 0, then one has p + 1 differential equations of order n' (and none of 
higher order). Thus, in this case all differential equations have the same order. 

Corollary 5.2. Let the differential equation (jS.ip be given where n = n'{p + 1) + n" with 
< n" < p + 1. Then (15. ip is equivalent to a set of n" differential equations of order n' + 1 as 
well as p — n" + 1 differential equations of order n' for the components fk{x). In particular, the 
set of solutions of (15. ip is a complex linear space of dimension n. 



Proof. It remains to check the dimension of the space of solutions. However, since the linear 
differential equation of order n' (resp. n' + 1) has n' (resp. n' + 1) linear independent solutions, 
one has in total n"{n' -\-l) + {p — n" + l)n' = n'[p+ 1) + n" = n linear independent solutions. ■ 

Let us check that the above equations given in Proposition 15.11 do indeed reduce to the ones 
given in (j5.4p - (l5.5p for the case n < p, i.e., n' = and n" = n. Consider the first equation 
of Proposition 15.11 Since n' = the right-hand side vanishes. It follows that co[n]q\fn{x) = 
— ^^Zq Cn~uW]qlfu{x)- Recalling cq = 1 and introducing i = n — this is equivalent to 

fn{x) = — ^ ^ Ci Y~i ''fn—i{x) 
u=0 ^^J'?- 
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which is the case = of (j5.4p . In the second equation of Proposition 15.11 the right-hand side 
also vanishes, implying "^^^QCn-ui^ + '^jqUl+ui^) = 0. Using cq = 1 and introducing i = n — v 
yields as above 

fl+n{x) =-2^Ci ^^^^j^i fl+n-iix) 

which are the (remaining) cases 1 < / < p — n of ()5.4p . In the third equation of Proposition 15. II 
only the first sum remains, implying 

n p— / 

^ Cn-u\l + 1^ -{p+ l)]q\dxfl+u^{p+l){x) = - ^ Cn~u[l + l']q\ fl+u{x) ■ 

Singling out again the summand z/ = n, recalling cq = 1 and introducing i = n — v yields 

/+n-(p+l) 

[l + n-{p+ l)]q\dxfi+n~{p+i){x) = - ^ Ci[l + n-{p+l)- i]q\dxfl+n~{p+l)~i{x) 

i=l 
n 

- X] Ci[l + n-i]qlfi+n-i{x) 
i=n+l~p 

which is exactly (IS.Sp . Thus, the system of equations given in Proposition 15.11 reduces for n < p 
indeed to (|5.4p - (j5.5p . Now, let us turn to the proof of the proposition. 

Proof of Proposition 15.11 The proof is straightforward but slightly tedious since one has to 
be very careful with the indices involved. Let us introduce cq := 1 as well as dg := c„_s for 
< s < n. Recalling furthermore that we can split n in a unique fashion as n = n' (p + 1) + n" 
with < n" < p, we can write the left-hand side of (jS.ip in the following form 

n n'-l p n" 

Inserting the expansion f{x, 6) = fo{x) + fi{x)9 + • • • + fp{x)6^ and using (|2.6p yields 

n n' — l p v~l r, 1 I 

s=0 f,=0u=0k=0 ^)+K\q. 

n'-l p p „ , I 



u=0 k=0 
n" p 



+ E E dn'{p+l)+u ^^^^^^^j^j ^x fk{x)e'' 



u=0 k=u 

We now consider the terms on the right-hand side of this equation separately. Let us begin with 
the second term. Introducing I = k — i^, we switch from the indices k, v to I, v and write this 
term equivalently as 

p p~l n'-l [I + 1^] ! 

EE E dti{p+i)+,y r,i ,'' d^fi+u{x)e^. 

1=0 u=0 /x=0 
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Let us turn to the fourth term. Here we also introduce I = k — u and switch from the indices k, u 
to I, u. Treating the range of v carefuhy, this term can equivalently be written as 

p min{p-l,n") _|_ j^] ] 

1=0 u=0 ^^J?- 

In the first term we introduce I = p + 1 — {u — k); it ranges from 1 to p. Switching from the 
indices k, v to I, v, this term is equivalent to 

A A '^sr^ [l + u- (p+ l)]g! a^^ . 

2^ 2^ Z^d^ip+i)+u (^xJl+u-ip+l)(X)(:> ■ 

1=1 u=p+l-l fj.=0 ^ 

Treating the third term in the same fashion (carefully treating the ranges) yields 

l=p+l-n" u=p+l-l 

Inspecting the resulting expressions for the four terms we see that upon combining them in 
powers of 9 there are three ranges to consider: 1) The range I = 0; here one has a contribution 
from the second and fourth term. 2) The range 1 < / < p — n"; here one has a contribution from 
the second and fourth term as well as the first term. 3) The range p — n" + 1 < I < p\ here all 
terms contribute. We can, therefore, write 

n p—n" p 

Y,dsV'f{x,e) = Eo{x) + Fi{x)9' + Yl Gi{x)e\ 

s=0 1=1 l=p+l-n" 

where we have abbreviated the corresponding functions Eq{x), Fi{x) and Gi{x) as follows: 

p n' — 1 n" 

Eoix) = YY1 d^<p+^)+'^W-^xfu{x) + X]'^"'{p+i)+z^M<?!'9"7r.(2;), 
u=o ^l=o 1^=0 

p-l n'-l [/ -I- 1 I mm{p-l,n") \1 , ] ] 

i.=OM=0 u=0 ^^^1- 

P n'-l \l + u] ! 

+ E E df,{p+i)+u^=^d^fi+^{x), 



v=p+l—l ^1=0 
n" 

Gi{x) = Fi{x)+ ^^n'(p+i)+.^=i^9:'+V/±.(^). 

u=p+l~l 



Thus, ()5.ip is equivalent to Eq{x) = 0, Fi{x) = = Gi{x) for all I. Let us start with Eq{x) = 0. 
This is equivalent to 

n" n' — l p 

E^"'{p+i)+'^M9!'9^7i/(2;) = - E E^M(p+i)+^M'?''5x/i^(a^)- 

u=0 /i=0 u=0 

Recalling ds = Cn-s as well as n = n'{p + 1) + n", this is equal to 

n" n' — l p 

Ycn"-uWW9x' fuix) = - Y^^{n'-f^){p+l)+{n"-i')WWdiiffix), 
u=Q 11=0 u=Q 
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which is the first asserted equation. Let us turn to Fi{x) = 0. The highest order of a derivative 
appearing is n'. Thus, switching also to the coefficients Cg, we can write Fi{x) = (with 
1 < I 1^ P — n") equivalently as 

n" n'-l / p-l 

^Cn"-u[l + iy]q\d^' fl+u{x) = " J] ^ ^{n' -,i)ip+l)+{n" -'^)['' + ^U^^ fl+iyi^) 
u=0 11=0 \iy=0 

P \ 

u=p-l+l / 

which is the asserted second equation. It remains to consider Gi{x) = 0. Here derivatives of 
order n' + 1 appear. Thus, switching also to the coefficients c^, we can write Gi{x) = (with 
p — n" + 1 <l <p) equivalently as 

n" p-l 

c„»_,[[ + i.],!a,"'+V/+.(x) = -Yl ^""-[^ + ^V-dxfi+u{x) 

u=p-[-l-l u=0 

n'-l / p-l 

- Yc(^n'-ti){p+l)+{n"-v)[l + l^Vd^fl+v{x) 

11=0 \ u=0 

P \ 



+ X] ^{n'-ti){p+l)+{n''-u)\L±J^W-d!ifl±u{x) 
u=p-l+l / 



which is the asserted third equation. 



6 Linear systems of first order 

Using the standard procedure to transform a linear differential equation into a system of 
differential equations of first order, we define v\{x^ 9) := /(x, 9) as well as Vj{x, 9) := T)^~^f{x, 9) 
for 2 < j < n. Forming the vector function u (x, 9) := {vi{x, 9), . . . , Vn{x, 9)Y, we can transform 
the differential equation ()5.ip into the form T>v{x,9) = Av{x,9) where the matrix A G C"'" is 
given in the usual form. Clearly, one is interested in the general system 

Vw{x,9) = Aw{x,9), (6.1) 

where the matrix A S C"'"' is arbitrary and w{x, 9) = {wi{x, 9), . . . , Wnix, 9)Y is a function to be 
determined. What one is in fact looking for is the generalization of the function eq{a^^^x; a9) = 
^aP+^x ^P^^ satisfying Veq{aP~^^x; a9) = aeq{a'P~^^x; a9) to the case with matrix arguments. 

Proposition 6.1. Let A G C"'" he given. For a given p > 1 define the matrix exponential 
function Eq by 



Eq{AP-^'x;A9) :=e^''^'^Yl 



{A9)'' 



,=0 t^]^' ■ 

The so defined function Eg has the following properties: 

1. In the scalar case n = 1 it reduces to the exponential function Cq. 

2. // the matrix A has diagonal form, i.e., A = diag(ai, . . . then one has 

Eg{AP+^x; A9) = diag {e,«+^x; ai9), eq{aP+^x; a^^)} . 
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3. It satisfies VEg{AP+^x; AO) = AEg{AP+^x; AO). 

4. It can also be written as Eg{AP+^x; A6) = YX=o ^Jy^e^''^'"'- 

Proof. The proof of the first property is clear. The second property is shown by a straightfor- 
ward computation and the third property follows as in the scalar case by a direct calculation 
from the definitions. The last property follows since powers of A commute. ■ 



Let us now consider the general linear system (j6.ip . Recalling the third property of Eq stated 
in Proposition 16.11 it is clear that w{x,9) := Eq{AP'^^x; A6)c - where c = (ci,...,c„)* with 



Cj € C is a vector of constants - is a solution of (j6.ip . It remains to be shown that no other 
solutions exist. This will be done by an explicit construction of the solution. 



Theorem 6.1. Let the linear system (16. ip with A £ C"'" be given. Any solution w{x, 9) of (16.1 
can be written as 



w 



{x,e) = Eq{AP+^x;A9)c, 



where c = (ci, . . . , c„)* with Ci G C is a vector of constants. In particular, the set of solutions is 
a complex linear space of dimension n. 

Proof. To show that all solutions can be written in this fashion we prove the asserted formula 
by constructing the solution. Inserting the function w{x, 9) = {wi{x, 9), . . . , Wn{x, 9)Y into (16. ip 
yields the following system 



1 " 
--—d^Wifiix) = y^auwi^pix), 
\p\g\ ^ 



1=1 



1=1 



which is equivalent to 



1 < i < n, 



1 < i < n. 



1 <k<p, 



f dxWifi{x) 


\ 




( wi^p{x) 


\ 












\ dxWn,o{x) 


/ 




\ Wn,p{x) 





1 wi^k{x) \ 


A 


' wx,k~i{x) \ 


\ Wn,kix) } 




V Wn,k-l{x) ) 


pie iteration shows that 


/ wx^k{x) \ 


j^k 


^ Wlfl{x) \ 


\ Wn,k{x) J 


1 " [k],\ 


\ Wn,o{x) ) 



l<k <p. 



l<k<p, 



and, consequently, 



/ dxWifi{x) 


\ 




( wi^o{x) \ 






= AP+^ 




\ dxWnfl{x) 


) 




\ Wnfl{x) ) 
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Thus, there exists a vector c = (ci, . . . , CnY of constants Cj € C such that 







( ^ 





















Collecting the above results shows that 
w{x, 9) = 




k=o 

as asserted. ■ 

Remark 6.1. It is interesting to note that there is no simplification of the proof or the resulting 
formula for the case p = 1 - in contrast to the scalar case where some simplifications occur. 

7 Some remarks concerning variable coefficients 

It is very important to note that we have here two types of variables, namely the ordinary 
"bosonic" one, i.e., x, and a nilpotent "fermionic" one, i.e., 9. The fermionic variable is more an 
algebraic object then an analytical one, but here the two variables are mixed. Let us consider 
the linear differential equation with variable coefficients, i.e., an equation of the form 

(P" + ci(x, + • • • + c„_i(x, 9)V + cn{x, 9)) fix, 9) = g{x, 9), 

where the coefficients Ci{c,9) are functions. We expect that it will be rather difficult to find 
solutions in general since the equation is more complex then the "ordinary" differential equation 
involving only the bosonic variable x. However, if we consider instead the coefficient functions to 
be depending only on the fermionic variable, i.e., Ci{x, 9) = Ci{9), then the resulting differential 
equation 

(P" + ci(^)P"-i + • ■ ■ + Cn-iie)V + c„(0))/(x, 9) = gix, 9) 

becomes more tractable since multiplication with the functions Ci{9) results only in a swirling 
of the coefficients - and no "ordinary" differential equations with variable coefficients have to 
be solved as subproblem. Instead of discussing these equations further we illustrate them with 
one of the most simple examples. 

Example 7.1. Let the differential equation 

{V-c{9))f{x,9) = Q (7.1) 

be given. Clearly, if c{9) were constant, i.e., c{9) = c, then we would have to solve the eigen- 
value equation 'Df{x,9) = cf{x,9) whose solution we have already found to be f{x,9) = 
Ceq{cP+^x; c9). Considering ^ and writing c(e)= ^f^^ ci9\ the ansatz /(x, 9)= Y,l=o fk{x)9^ 
yields the system 

I P ^ 

YY-^dxhix) = '^Cp-ifi{x), [k + l\gfk+i{x) = '^Ck-ifi{x), 0<k<p-l. 

^^1- 1=0 1=0 
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The strategy to solve this equation is clear: Iteration of the second equation yields a relation 
of the form fk{x) = Ck{c{6),p)fQ{x) which can be inserted into the first equation to obtain 
a differential equation dxfo{x) = D{c{9),p)fo{x) which can be solved as fo{x) = Ce^^^^^^'^^^ . 
Inserting this into the result for fk{x) yields then fk{x) and, consequently, f{x,9). However, to 
determine the concrete values for the coefficients Ck{c{6),p) and D{c{9),p) seems to be difficult 
in the general case. As an illustration we consider the special case p = 2. Here one obtains for 
k = the relation fi{x) = cofo{x) and for k = 1 the relation [2]g/2(x) = ci/o(x) + 00/1(2;), 

implying f2{x) = fo{x)- These results can be inserted into the first equation and give 
d.foix) = {cl + (1 + [2],)coCi + [2],C2} /o(x) = D{c{e), 2)/o(x). 

Thus, foix) = Ce^^^^^^''^^^ . Inserting this into the equations for the components fkix), we have 
finally found 

f{x,e) = Ce^W^)-2)- |l + co^ + ^ll^^'} • 

Note that in the case where c{9) is constant, i.e., c{6) = c, one has D{c,2) = = c^+^ and 
the solution given above simplifies to f{x,0) = Ce'^^^^^il + c9 + ^^V) = C eq{c^~^^ x] c9) , as it 
should. 



Remark 7.1. As noted already in the introduction, a certain differential equation involving a 
para-Grassmann variable was solved in [58] (by writing f{x,9) = X^^=q /fc(a;)^'^, inserting this 
into the differential equation and comparing coefficients of 9^). Adapting the notation to the 
one used in the present paper, the differential equation considered in [38] can be written as 



P (_n\l 

9. + (^x + ^) ^ i-^ai 



/=0 



f{x,9) = \f{x, 



where A, G C. Clearly, the structure of this differential equation differs from the structure of 
the differential equations considered in the present paper - see, e.g., equation p.4p for a simple 
eigenvalue problem in our context. 



8 A simple family of nonlinear differential equations 

In the above sections we have considered certain types of linear differential equations and 
have found them to be very similar to "ordinary" differential equations (with only "bosonic" 
variables). In this section we want to discuss some equations which are closer to partial 
differential equations in "bosonic" variables. Let us start with the very simple nonlinear equa- 
tion (^/(x, 9))'^ = 0. Solving this equation is straightforward: All one has to do is to insert (j2.3p 
and compare coefficients. The result is that the general solution is given in the case p = 1, 2 by 
fix, 9) = fo{x) and in the general case of p > 3 by f{x,9) = fo{x) + X]^=2+[h] f'^^^)^''' where 
we have denoted by [x] the largest integer smaller or equal than x. Thus, the general solution 
depends on p — [|] functions. More generally, one may consider (P"'/(x, 0))" = for arbitrary 
n G N and m < p. Writing V"^f{x,9) = Ylf=o^li^)^^ ^^'^ considering (X]]Lo ^'(^)^')" ~ ^ 
implies the vanishing of ho{x), . . . , /ijpj (x). Applying (|2.4|) shows that hQ{x) = fm{x) and that 

the other hj{x) are given in a similar fashion where an interesting connection between m, n 
and p appears. These results are collected in the following theorem. 



On Linear Differential Equations Involving a Para-Grassmann Variable 



23 



Theorem 8.1. For arbitrary p, n G N and m < p the general solution of the differential equation 
(V^ f , 9))'"' = is given in the case 1 < m < p — [^] by 

m—l p 

f{x,e) = Y,fk{x)e^+ E /'^•(^)^' 

and in the case p — < m <p hy 

fc=0 fc=m_(p_[Jl]) 

where the appearing functions fk{x) are arbitrary [note that there are in both cases p — [^] of 
them) and the Cfc € C are constants. Thus, the set of solutions has infinite dimensions. 

Using the "periodicity" described in (12. 5p . the case where m > p + 1 can be reduced to the 
one considered in the above theorem. In fact, writing m = m'{p + 1) + m" with < m" < p, 
one has V"'f{x,d) = V"'"h{x,9) with h{x,e) = f{x,e). 



9 Conclusions 

In this paper we have considered linear differential equations involving a para-Grassmann vari- 
able vanishing in {p + l)-th order. Suppose that we are given an equation of order n = 
n'{p + 1) + n" . Expressing the parasuperfunction f{x,6) through the components fk{x) and 
inserting this into the differential equation results in a system of n" differential equations of 
order n' + 1 as well as p — n" + 1 differential equations of order n' for the components fk{x). 
In the case n < p one has, therefore, a system of p — n + 1 recursion relations and n differential 
equations of first order. The key to the solution of this system consists in expressing the 
"higher components" through the "lowest components" and then solving the resulting system 
of ordinary differential equations for the "lowest components" . Having found the explicit form 
of the "lowest components" then allows one to describe the "higher components" and, conse- 
quently, f{x,9) itself. For the linear differential equation with constant coefficients of order 
n < p we have determined its solutions and found that the set of solutions is a complex lin- 
ear space of dimension n (the latter result also holds for arbitrary n). Characteristic for these 
linear differential equations with constant coefficients of order n < p is an interesting blend of 
"ordinary" differential equations (involving only the "bosonic" variable x) and recursion rela- 
tions. These recursion relations are intimately connected to the Fibonacci numbers and their 
generalizations to higher order. It is interesting to note that in the results there is no difference 
between the cases p = 1 (i.e., an "ordinary" Grassmann variable) and p > 1, although in the 
calculations (and proofs) there are striking differences: In the case p = 1 the above-mentioned 
recursion relations are empty (i.e., not existent) and only a system of differential equations has 
to be solved. It also deserves to be mentioned that the appropriate exponential function is of 
utmost importance for describing the solutions. This is similar to the "bosonic" case but there 
exists a decisive difference: Here the "derivative" V is not a derivation, i.e., one does not have 
the usual law for the derivative of a product at hand. In particular, in the case corresponding 
to degenerated eigenvalues it is a priori not clear what the analogue of the usual ansatz x^e^^ 
should be in the para-Grassmann case and one of the main difficulties met in the present work 
was to solve also these cases (even with hindsight - i.e., the formulas for the solution at hand - 
it is not clear whether one would have made such an ansatz). Following the usual approach, 
the linear differential equation with constant coefficients can be transformed into a system 
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of differential equations of first order. For such a general system we determined its solution 
and showed in particular that the set of solutions is also a complex linear space (of expected 
dimension) . 

We also considered some other classes of differential equations which, however, seem not 
as well behaved as the classes considered above and, therefore, deserve a closer study in the 
future. As a first example we mention the scalar linear differential equations where variable 
coefficients are allowed. We considered only an example of the case where the coefficients are 
functions of the para-Grassmann variable 9 (no dependency on x). Again, one finds a blend 
of "ordinary" differential equations and recursion relations but due to the variable coefficients 
the components in the recursion relations obtain a "twist", resulting in slightly cumbersome 
expressions. It is obvious that for the general case some new ideas are needed. As a second 
example we considered a very simple family of nonlinear differential equations. Due to the 
nilpotency of the para-Grassmann variable 9 we found a behavior than differs drastically from the 
linear case: Although the ususal blend of recursion relations and ordinary differential equations 
appears there exist (for appropriate choices of the involved parameters) several components of 
the solution for which no restrictions have to be satisfied. Thus, the dimension of the set of 
solutions is in general not m (as one would have expected for the differential equation of order m) 
but infinite! Thus, the close analogy to the ordinary ("bosonic") differential equations found in 
the linear cases breaks down and one finds a behaviour which reminds one more of "bosonic" 
partial differential equations. Clearly, this is also a point where a closer study might reveal 
further interesting connections. 

Before closing let us remark that - in addition to the above-mentioned points concerning 
differential equations involving one "bosonic" variable x and one para-Grassmann variable — 
one can also consider partial differential equations involving several bosonic and para-Grassmann 
variables. Since in this case the para-Grassmann variables are ~ in addition to being nilpotent - 
noncommutative, it is to be expected that it will be much harder to obtain results for these para- 
Grassmann partial differential equations than for the "ordinary" para-Grassmann differential 
equations considered in the present paper. Before beginning this journey into the unknown 
territory full of treasures one should prepare oneself to meet surprises - pleasant or otherwise - 
around every corner and it seems mandatory to become comfortable with those of them showing 
up already in the one-variable case. This is what we began in the present paper. 
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